Research in Astronomy and Astrophysics, 24:025016 (Spp), 2024 February 
© 2024. National Astronomical Observatories, CAS and IOP Publishing Ltd. Printed in China and the U.K. 


CIV Broad Absorption Lines with Different Ionization Levels 


Ying-Ru Lin’, Cai-Juan Pan’, and Wei-Jian Lu! © 
' School of Information Engineering, Baise University, Baise 533000, China; yingru_lin@qq.com, william_lo@qq.com 
School of Materials Science and Engineering, Baise University, Baise 533000, China 
Received 2023 July 15; revised 2023 September 19; accepted 2023 October 17; published 2024 January 31 


Abstract 


This paper explores the correlation between the fractional variation of the ionizing continuum and CIV broad 
absorption lines (BALs) with different ionization levels. Our results reveal anti-correlations between fractional 
variation of the continuum and fractional equivalency width (EW) variation of the C IV BALs without Al 1u BAL/ 
mini-BALs at corresponding velocities, providing evidence for the widespread influence of the ionizing continuum 
variability on the variation of HiBALs. Conversely, for CIV BALs accompanied by Al It BAL/mini-BALs 
(LoBAL groups), no significant correction is detected. The absence of such a correlation does not rule out the 
possibility that variations in these low-ionization lines are caused by ionizing continuum variability, but rather 
suggests the influence of BAL saturation to some extent. This saturation effect is reflected in the distribution of the 
fractional EW variation, where the CIV BAL group accompanied by Al m BAL has a smaller standard deviation 
for the best-fitting Gaussian component than the two BAL groups without Al m BAL. However, the distribution of 
fractional variation of their continuum does not show any significant difference. Besides the saturation influence, 
another potential explanation for the lack of correlations in the LoBAL groups may be the effects of other 
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Relations between the Fractional Variation of the Ionizing Continuum and 


variability mechanisms besides the ionization change, such as clouds transiting across the line of sight. 
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1. Introduction 


Quasar outflows have been proposed as a kind of potential 
active galactic nucleus (AGN) feedback mechanism, which 
could play a role in terminating star formation in the interstellar 
medium (ISM) and regulating accretion onto the supermassive 
black hole (SMBH; Fabian 2012; King & Pounds 2015). In 
rest-frame UV spectra of quasars, outflows typically appear as 
broad absorption lines (BALs), with absorption widths greater 
than 2000 km s~! (Weymann et al. 1991). “Mini-BALs” are 
defined as absorption troughs with full widths at half minimum 
from a few hundred to 2000 km s~! (Hamann & Sabra 2004). 
Most BALs have only high-ionization transitions, such as N V, 
SiIV, and CIV (referred to as HiBALs, as mentioned by 
Weymann et al. 1991). A smaller fraction, about 15% in BAL 
quasars, also exhibits absorption from low-ionization transi- 
tions, including AlN, All and Mg II (known as LoBALs, as 
reported by Voit et al. 1993 and Gibson et al. 2009). Another 
subgroup, known as FeLoBALs, represents a small subset of 
LoBALs that exhibit absorption of Fell or Fem (Wampler 
et al. 1995). 

BAL variability can serve as a means to assess BAL 
structures, wind locations, and dynamics, and constrain the 
physical mechanisms responsible for AGN outflows. Time 
variability in individual sources or samples with multi-epoch 
observations has been reported for HiBALs (e.g., Capellupo 


et al. 2011, 2012; Filiz Ak et al. 2013; Yi et al. 2019a, 2019b; 
He et al. 2019), for LoBALs (Zhang et al. 2011; Vivek et al. 
2014; McGraw et al. 2017; Sun et al. 2017; Yi et al. 2019a; Yi 
& Timlin 2021), and for FeLoBALs (Vivek et al. 2012; 
McGraw et al. 2015; Zhang et al. 2015; Rafiee et al. 2016; 
Stern et al. 2017). The observed BAL variability can usually be 
explained by variations in the ionization parameter (e.g., 
Barlow 1994; Crenshaw et al. 2003; Filiz Ak et al. 2013; Wang 
et al. 2015; He et al. 2017; Lu & Lin 2018; Lu et al. 2018; 
Vivek 2019), clouds transiting across the line of sight (e.g., 
Lundgren et al. 2007; Gibson et al. 2008; Hamann et al. 2008; 
Krongold et al. 2010; Hall et al. 2011; Filiz Ak et al. 2012; 
Vivek et al. 2012; Capellupo et al. 2014; Shi et al. 2016; Vivek 
et al. 2016; Capellupo et al. 2017), or a combination of them 
(e.g., Capellupo et al. 2012; Vivek et al. 2014). One effective 
way to determine the variation mechanism of BAL in quasar 
sample is to investigate the correlation between variability in 
absorption lines and that of the continuum. Early study based 
on relatively small quasar samples reported no detectable 
correlation between changes in the BAL equivalency width 
(EW) and the continuum flux (Gibson et al. 2008; Vivek et al. 
2014; Wildy et al. 2014). However, recent studies have found 
anti-correlations between the fractional variation of the 
continuum and that of both CIV and SiIvV BAL troughs in 
several BAL quasar samples (Lu & Lin 2018; Lu et al. 2018; 
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Vivek 2019). These anti-correlations reveal the widespread 
influence of ionizing continuum variability on BALs, providing 
evidence for photoionization-driven BAL variability. 

The question of whether these correlations extend to lower 
and higher ionization levels is still under debate. Vivek et al. 
(2014) found no clear correlation between absorption line and 
the continuum flux variabilities from a quasar sample of 22 
low-ionization quasars that include MgII and AIM BALs. 
However, Vivek et al. (2014) pointed out that targets with large 
absorption line variability also present large variability in their 
light curves, and they concluded that the observed variability 
can be well explained by the combination of two mechanisms: 
one is the continuum variations, and another is the clouds 
transiting across our line of sight. 

In typical disk-wind models (e.g., Murray et al. 1995; Proga 
et al. 2000; Higginbottom et al. 2013), BALs with different 
ionization levels can be explained by viewing angle effect. 
More concretely, lower degrees of ionization are observed 
along lines of sight closer to the accretion disk plane (Hamann 
& Sabra 2004; Baskin et al. 2013, and Figure 16 of Filiz Ak 
et al. 2014). This view is supported by observations as well as 
simulation studies. For example, Filiz Ak et al. (2014) studied 
the variability characteristics of CIV BAL troughs and whether 
they also contained accompanying SiIV and/or Alm BAL 
troughs. The measurements that they conducted demonstrate 
the correlation between the changes in ionization level, 
kinematics, and column density, which is consistent with the 
prediction of accretion-disk wind models. For example, low 
ionization absorption features in line of sight typically exhibit 
deeper and broader CIV BAL troughs, but exhibit smaller 
minimum velocities and less variation. Recent numerical 
simulations have also shown that when low ionization 
absorption features presented, the BALs are wider and deeper, 
whereas high ionization species exhibit higher blue-edge 
velocities compared to the low ionization lines (Matthews 
et al. 2016). 

Filiz Ak et al. (2014) have reported the correlated changes of 
BAL with different ionization levels. Building upon their work, 
this paper aims to explore the correlation between the fractional 
variation of the ionizing continuum and CIV BALs with 
different ionization levels. The rest of this paper is organized as 
follows. Section 2 presents the data preparation and measure- 
ments of BALs and spectra. Section 3 contains the results, 
discussions, and conclusions. 


2. Data Preparation and Measurements 


The initial sample of this paper containing 671 quasars with 
at least two-epoch observations was selected by Filiz Ak et al. 
(2014) from SDSS-I/II/IIL. The selection criteria for these 671 
targets were described in Section 2 of Filiz Ak et al. (2014). In 
their study, Filiz Ak et al. (2014) identified CIV, SiIV and 
Alt BAL/mini-BALs in the spectra of these quasars and 
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Table 1 
Number and Spearman Rank Correlation Analysis Results of Nine Groups 
of BALs 
Total BAL Variation BAL 

BAL types Number* Number” P Ps 
C IVoo 113 62 —0.46 1.97E-04 
CIV 271 161 —0.43 1.25E-08 
C IVso 246 144 —0.26 1.33E-03 
Si IVso 246 135 —0.39  2.67E-06 
C IVsa 32 16 —0.29 0.27 
C IVsa 94 49 —0.26 0.07 
C IVsa 95 38 —0.25 0.13 
Si IVsa 95 39 —0.18 0.27 
Al Msa 95 32 —0.31 0.08 
Notes. 


è BAL number selected by Filiz Ak et al. (2014). 

BAL number selected in this article. 
4 The Spearman coefficient of correlation (r) and P-value (P) are given in 
bold where P <0.01. For the C IVoo, C IVso, C IVso, C IVsas C IVsa, and C IVsa 
groups, the AEW/ (EW) is measured for C Iv BAL troughs. For the Si IVso 
and SiIVgq groups, the AEW/ (EW) is measured for Si Iv BAL troughs. For 
the Al Msa group, the AEW/ (EW) is measured for Al M BAL troughs. 


classified the CIV BALs into six groups according to whether 
there are SiIV and/or Alli! BALs/mini-BALs at corresp- 
onding BAL regions. Specifically, “CIVo9” represents CIV 
BALs without accompanying SiIV and Alm BALs at 
corresponding velocities during both epochs; “C IVsoọ” repre- 
sents CIV BALs with a Si 1V BAL at corresponding velocities 
in one or both epochs, but no corresponding Alm BAL; 
“C IVsa” represents CIV BALs with both accompanying Si 1V 
and Al m BALs in either epoch; “C IV,9” represents CIV BALs 
with a Si IV mini-BAL, but no accompanying Al m BAL/mini- 
BAL is detected; “CIV,,” represents CIV BALs with both 
accompanying Si IV and Al M mini-BALs, “C IVs,” represents 
CIV BALs with a Si IV BAL as well as an Al Il mini-BAL. In 
addition, Filiz Ak et al. (2014) defined the corresponding Si IV 
BALs of C IVso as “Si IVso,” the corresponding SiIV and Al M 
BALs of CIVgaq as “SilVg,q” and “Al Msa,” respectively. It is 
worth noting that among these nine groups, CIV, 
CIV.9,C IVsq and SiIVso belong to HiBALs, while CIV,,, 
C IVsa, CIVgq, SiIVsa and Al Msa belong to loBALs. 

We calculated the EW variation (AW) of all these detected 
BAL/mini-BALs, adopting the EW values measured by Filiz 
et al. (2014) and applied Equation (1) from Lu et al. (2018). We 
then selected BALs/mini-BALs based on the criterion of 
AW > 3ø. The number of BALs included in each of the nine 
groups from Filiz Ak et al. (2014) and those selected for this 
article are listed in Table 1. We then calculated the fractional 
EW variation (AEW/ (EW)) of our selected CIV, SiIv, and 
Alm BALs using Equation (2) in Lu et al. (2018). 

To obtain the flux density values, we fitted a power-law 
continuum for each spectrum utilizing the same procedure in 
Lu et al. (2018). We then used the power-law continuum flux at 
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Figure 1. Fractional variation of BALs (AEW/ (EW)) as a function of the fractional variation of the continuum (AF eont/ (Fcont)) for nine groups of BALs. For the 
C IVoo, C IVg0, C IVgo, C IVs, C IVgq, and C IVs groups, the AEW/ (EW) is measured for C Iv BAL troughs. For the Si IVso and Si IVgq groups, the AEW/ (EW) is 
measured for Si IV BAL troughs. For the Al Msa group, the AEW/ (EW) is measured for Al M BAL troughs. 


1350 A from the two-epoch spectra of the same quasar to 
calculate the fractional variation of the ionizing continuum 
(AF cont/ (Feont)), following Equation (4) in Lu et al. (2018). 
Figure 1 depicts the fractional variation of BALs 
(AEW/ (EW)) against the fractional variation of the con- 
tinuum (AF .ont/ (Feont)) for each group, and the corresponding 
results of the Spearman rank correlation analysis are listed in 
Table 1. 


3. Results and Discussions 


As shown in Table 1, anti-correlations between the fractional 
flux variations of the continuum and fractional EW variations 
of the BALs in C IVoo, C IVso, C IVso and Si IVso groups have 
been confirmed by the Spearman rank correlation analysis 
(with P-values of Spearman coefficient less than 0.01). It is 
worth noting that the CIV BALs with lower ionization exhibit 
smaller correlation coefficients. When the CIV BALs are 
accompanying with Alm BAL/mini-BALs (CIV.a, CIVsa, 
CIVsa SiIVsa and AlMsa), no significant correction is 
detected. These statistical results are consistent with the 
expectations of the disk-wind model and provide important 
insights into the variation mechanism of BALs with different 
ionization levels. 


On one hand, we confirm the anti-correlations between the 
fractional variations of the continuum and BALs in four groups 
(CIVo9, CIVso, CIVso and SiIVso) of HiBALs. Previous 
studies have observed anti-correlations between the fractional 
variations of the continuum and absorption lines in several 
BAL samples (Lu & Lin 2018; Lu et al. 2018; Vivek 2019) as 
evidence for ionization-driven BAL variation. However, these 
studies did not compare the impact of different ionization states 
of BAL on the correlation. Our study expands on this by 
finding anti-correlations between the fractional flux variations 
of the continuum and fractional EW variations for four BAL 
groups without Alm BAL/mini-BAL, covering a wider 
ionization-potential range. This finding provides further 
evidence of the widespread impact of ionizing continuum 
variability on the variation of HiBALs (e.g., Weymann et al. 
1991). Additionally, in accordance with the typical disk-wind 
model (e.g., Murray et al. 1995; Proga et al. 2000; 
Higginbottom et al. 2013), if different groups of BALs 
represent different viewing inclinations, the anti-correlations 
observed in our study provide evidence for the ubiquitous 
effects of ionizing continuum variability across a wide range of 
viewing inclinations. 

On the other hand, we found no significant correlations 
between the fractional flux variations of the continuum and 
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Figure 2. AF.ont/ (Feont) and AEW/ (EW) distributions for the BALs in our CIV 9 (dotted—dashed blue), C IVso (dashed green), and C IVs; (solid red) groups. 


fractional EW variations for the five BAL groups with Al m 
BAL/mini-BAL (C IV.a, C IVsa, CIVga, SiIVsa and Al Iga), 
with P-values of Spearman coefficient greater than 0.01. We 
speculate that the lack of correlations in these groups may be 
due to the effects of BAL saturation or other variability 
mechanisms besides the ionization change, such as clouds 
transiting across the line of sight. 

The observations of BAL profiles suffering from saturation 
have been reported in several previous studies (e.g., Arav 1997; 
Hamann 1998; Arav et al. 2001; Leighly et al. 2009; Borguet 
et al. 2012; Xu et al. 2018). For instance, Filiz Ak et al. (2014) 
have shown that when BAL troughs from lower ionization 
transitions are present, CIV troughs tend to be stronger and 
wider (see their Figure 5) and exhibit more saturation, but less 
fractional EW variation (see their Figure 11). This suggests that 
CIV troughs might experience more saturation when BAL 
troughs from lower ionization transitions are present. Hamann 
et al. (2019) also obtained empirical results showing that 
compared to the HiBALs, the CIV BALs in the composite 
spectrum of the LoBALs are deeper and wider (see their 
Figure 5). Using two sub-samples differing on the absorption 
trough depth, Vivek (2019) found that the shallow trough sub- 
sample exhibits an even stronger correlation while the deep 
trough sub-sample appears to be no correlation between the 
BAL variability and the continuum variability (see their 
Table 3). These results point out that the saturation of the 
BALs has a considerable effect in the correlation between the 
BAL variability and the continuum variability. The saturation 
can also be inferred from our Figure 2, which shows that the 
AEW/ (EW) distribution of the CIVsq group has a smaller 
standard deviation (ø) for the best-fitting Gaussian component 
(o =0.213) compared to the CIVgg (o = 0.390) and C IVso 
(o = 0.357) groups. This is because for saturated troughs, the 
fractional EW variation measurements can only reflect the 
lower limits of the true variations in optical depth and column 
density, thus weaken the correlation between the absorption 
line variability and continuum variability. However, the 
AF cont/ (Feon) distribution of the CIVs,q does not show any 
significant difference compared to the C IVoo and C IVgo groups 


(according to the K-S test results, with P > 0.026). Considering 
the significant impact of line-saturation in the CIVgq group, we 
speculate that the effects of ionizing continuum variability may 
also exit in the BALs that accompany Al m BAL/mini-BALs. 
In other words, the ubiquitous effects of ionizing continuum 
variability may extend to even lower viewing inclinations. 

Another potential explanation for the lack of correlations in 
our LoBAL groups could be the transiting of the clouds across 
our line of sight. Transiting scenario has been reported to be a 
feasible cause of the absorption-line variability (e.g., Hamann 
et al. 2008; Capellupo et al. 2017). For example, if a highly 
saturated CIV BAL exhibits clear variability, it provides 
evidence in favor of the transiting scenario. This is because a 
highly saturated CIV BAL is not easily affected by small 
changes in ionization, but can easily exhibit variations due to 
the transiting of outflows (e.g., McGraw et al. 2018). The 
presence of a PV BAL can serve as a detector for the highly 
saturated CIV BAL in the corresponding velocity (e.g., 
Capellupo et al. 2014, 2017; McGraw et al. 2018). To put 
further constraints on the origin of variability of our LOBAL 
subsamples, it is necessary to check individual spectra, which is 
beyond the scope of this paper. 

In summary, we have discovered anti-correlations between 
the fractional variation of the continuum and fractional EW 
variation of the CIV BALs in four HiBAL groups, while no 
significant correlation was found in five LoBAL groups. The 
anti-correlations presented in four HiBAL groups demonstrates 
the widespread impact of ionizing continuum variability on the 
variation of outflow absorption. The lack of correlations in five 
LoBAL groups does not rule out the possibility that variations 
in these low-ionization lines are caused by ionizing continuum 
variability, but rather suggests the influence of BAL saturation 
to some extent. We have tested the AEW/ (EW) and 
AF cont/ (Fcon) distributions of the C IVoo, C IVso, and C IVsa 
groups (Figure 2) to reveal the saturation influence. Another 
potential explanation for the lack of correlations in the LoBAL 
groups may be the effects of other variability mechanisms 
besides the ionization change, such as clouds transiting across 
the line of sight. 
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